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Abstract The mathematical theory of gravitational lensing has revealed many generic and global prop- 
erties. Beginning with multiple imaging, we review Morse-theoretic image counting formulas and lower 
bound results, and complex-algebraic upper bounds in the case of single and multiple lens planes. We 
discuss recent advances in the mathematics of stochastic lensing, discussing a general formula for the 
global expected number of minimum lensed images as well as asymptotic formulas for the probability 
densities of the microlensing random time delay functions, random lensing maps, and random shear, and 
an asymptotic expression for the global expected number of micro-minima. Multiple imaging in optical 
geometry and a spacetime setting are treated. We review global magnification relation results for model- 
dependent scenarios and cover recent developments on universal local magnification relations for higher 
order caustics. 

Keywords Gravitational lensing • singularities 



1 Introduction 

1.1 Overview and Conventions 

Two important anniversaries related to gravitational lensing occurred in 2009: ninety years ago the 
first observation of this effect was announced at a joint meeting of the Royal Society and the Royal 
Astronomical Society, as a successful test of Einstein's new theory of gravity; and thirty years ago 
Walsh, Carswell and Weyman reported the first observation of an extragalactic example of lensing. 
Especially since then, the subject has become a thriving research field at the interface of astronomy, 
theoretical physics and mathematics. Some current research highlights on astrophysical and cosmological 
applications of lensing have been discussed earlier in this Special Issue, as well as possible lensing tests 
of modified theories of gravity in the spirit of the original corroboration of General Relativity. 

Of course, it has also emerged that gravitational lensing theory is a rich research area in its own right 
within mathematical physics. This aspect can be approached from three different directions: the widely 
used and astrophysically important thin-lens, weak-deflection approximation^ optical geometry, which 
considers the properties of spatial light rays, a simplification that also makes the method applicable 
to astrophysically relevant models; and a full general relativistic spacetime method that studies null 
geodesies. These approaches have proved to be mathematically quite rich, with applications of singularity 
theory, differential topology, Lorentzian geometry, algebraic geometry, and probability theory. In this 
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Fig. 1 A schematic of single-plane gravitational lensing. A pointlike light source is at y on the light source plane S. A 
light ray from the source is deflected through an angle a by the gravitational influence of the lens on the lens plane L. 

review article, we discuss recent work in this direction on two aspects of the weak deflection limit, image 
counting and magnification. 

One of the most basic problems in gravitational lensing is the number of images produced. Yet, already 
this apparently simple question turns out to be difficult. Image counting results using Morse theory and 
complex methods are reviewed in Section[3]for the single lens plane case, and in Section[3]for multiple lens 
planes. The expected number of images in stochastic lensing is discussed in Section^ in particular for the 
asymptotic microlensing case of Section [5j From the point of view of optical geometry, image multiplicity 
is also a global effect as outlined in Section [6l Finally, conditions for the occurrence of multiple images in 
spacetime are summarized in Section[7l Going beyond image number, finer information on lensed images 
can be gained from the magnification. It turns out that image magnifications obey global magnification 
relations for certain lens models, and local magnification relations near singularities up to and beyond 
codimension three, which are universal. These results, which involve deep properties of singularities and 
algebraic geometry, are discussed in Sections O [Qj and [TOl respectively. 

Theorem Convention: Our criteria for deciding whether a mathematical result is called a theorem will be 
driven by its importance in advancing our understanding of the physical and/or mathematical aspects 
of lensing. 

Citation Conventions: The first time a result from a paper is mentioned, the name(s) of the author(s) are 
stated next to the bibliographic reference. Thereafter, citations of the same result have the bibliographic 
reference without the names. The authors of this article are cited as AOP for Fetters and MCW for 
Werner. 

Miscellaneous Conventions: All light sources are treated as point-like. The symbol = indicates that a 
definition is given. Physical quantities are given in appropriate dimensionless forms. 

1.2 Notation for Some Basic Weak Deflection Lensing Concepts 

Denote the dimensionless potential of a gravitational lens by ip- The (dimensionless) surface mass density 
K and magnitude F of shear are: 



where x = (u, v). Note that x is dimensionless, i.e., x = r/di, with r the physical impact vector in the 
plane of the physical lens at angular diameter distance d^; see Figure [TJ Since ^ |r|, we can treat x 
as an angular vector. A singularity of -0 is a point a G such that either —oo or V^?A(x) — ^ oo 

as x ^ a. In particular, an infinite singularity of V' is a singularity a for which ?A(x) ~^ —oo as x — > a. 
We shall call tp nonsingular if it has no singularities. The potential ip is assumed to be smooth (C°°) 
everywhere on M^, except on the set A of singularities of i/j. 

Let Ty : L — > R be the (dimensionless) single-plane time delay function induced by a lens potential 
ip, where L — M."^ — A with A the set of singularities of V': 



Unless stated to the contrary, we assume that yl is a finite set. The point y lies in the light source plane 
5' = ; see Figure [TJ The light source plane S is the set of all "angular" source positions y = s/ds, 



«:(x) = -VV(x), r'(x) = -(V^uu(x) - Vvv(x))' -t- 02^(x) 
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where s is a physical (hnear) source position in the Euchdean plane at angular diameter distance ds and 
orthogonal to the line of sight. 

By Fermat's principle [721164] . the light rays connecting y to the observer are given by the critical 
points of Ty, i.e., solutions x G i of 

= VTy(x) = -y + x- V^(x), 

where the gradients are with respect to the rectangular coordinates x = (u, v) . The above equation 
determines the single-plane lensing map corresponding to Ty : L — > K, namely, the transformation 
rj : L ^ S defined by: 

J7(x) = x — V?/'(x). 

The lensed images of a light source at y are the elements of J7^^(y) or, equivalently, the critical points 
of Ty. Through this correspondence, we speak of minimum, saddle, and maximum lensed images. 

Example (Microlensing) : We define microlensing generally as due to a lens consisting of g stars with 
masses mi, . . . , at respective positions ^i, . . . , ^g, continuous matter with constant density Kc > 0, 
and an external shear 7 > 0. The lens potential of microlensing is then given by: 

V.,(x) = ^|xp - ^(u^ - v2) + ^ m, log |x 
where x — (u, v). The induced time delay function Tg y at y is 

T,,y(x) = i|x - yp - ||xp + ^(u^ - v^) - E l°g 1^ - ^.1 



and the corresponding lensing map r] is 
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^sW = ((1 - + 7)u, (1 - Kc - 7)v) - E 



x-|. 



We shall refer to the lensed images in microlensing as micro-images and may even speak of micro-minima 
to designate minimum images in this context. 

The {absolute) magnification of a lensed image is physically the ratio of the flux of the image to the 
flux of the light source, which is given mathematically by (e.g., [64,, P- 85): 

Mag(x; y) = , , , J7(x) = y. 

I det[Jac i7j(x)| 

Note that det[Jac 77] (x) = det[Hessry](x) = (l — k(x))^ — r'^{-x). The signed magnification of lensed 
image x^ of y is defined by: 



(-l)ind(x.)Mag(x,;y), r7(x,) 



Here Ind(xi) is the Morse index of x^, i.e., the number of negative eigenvalues of (HessTy)(xi), which is 
0, 1, and 2 for a minimum, saddle, and maximum, respectively. 

The set Crit(j7) of critical points of t/ is the set of all (formally) infinitely magnified lensed images 
for all light source positions in S, which consists of all x € L where det[Jac T7](x) = 0. When Crit{ri) 
consists of curves, we shall speak of the critical curves of rj. The set of caustics of rj is Caustic(j7) = 
?7(Crit(?7)) C S, which is the set of all light source positions from which there is at least one infinitely 
magnified lensed image. For physically relevant settings, the set Crit(rj) is bounded, which yields that 
Crit(?7) and Caustic(j7) are compact ([M], P- 293). 

The single-plane time delay function Ty is said to be subcritical at infinity if for each non-caustic 
point y and |x| sufficiently large, the eigenvalues Ai(x;y) and A2(x;y) of HessTy(x) are positive. This 
means that 

Ai(x;y) = l-«:(x)+r(x)>0, Ai(x; y) = 1 - «:(x) - ^(x) > 0. (1) 
Condition ([T]) implies: 

< k(x) < 1, < r(x) < 1. 
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In addition, the time delay surface (graph) of Ty then has positive Gauss curvature for |x| sufficiently 
large. When Ty is subcritical at infinity and ry(x) ^ oo as |x| oo, we call Ty isolated. 

Given that the lensed images of a light source at y are in 1-1 correspondence with the critical points 
of Ty, we shall characterize each lensed image as either a minimum, saddle, or maximum. 

Notation: 

1. N — total number of images. 

2. Nniin — number of minimum images. 

3. iVsad = number of saddle images. 

4. iVmax = number of maximum images. 

5. N+ — Nynin + ^max = number of positive parity images. 



2 Multiple Imaging in Single-Plane Lensing 

2.1 Image Counting Formulas and Lower Bounds: Single Plane Lensing 

Einstein determined in 1912 [7T] that a lens consisting of a single star will produce two images of a 
background star that is not on a caustic. Two stars on the same lens plane can produce either three or 
five images for light sources off a caustic, a result found in 1986 by Schneider and Weiss [7^. The method 
in [73], however, involved lengthy calculations that directly solved the lens equation for the images, an 
approach that would be impossible to carry over to any ffnite number of stars or more general mass 
distributions. 

In 1991, AOP ^61] approached the image counting problem by employing Morse theory under bound- 
ary conditions and generic properties of time delay functions to obtain a general theorem yielding counting 
formulas and lower bounds for the number of images. The Morse theoretic approach is particularly pow- 
erful because it gives specific counting information about the number of images of different types and 
extends naturally to fc-lens planes and a general spacetime setting. For simplicity, we state the results in 
[5T] only for generic subcritical lensing (see [BH Chap. 11] for more): 

Theorem 1 [61] (Single-Plane Image Counting Formulas and Lower Bounds) Let Ty : L ^ M. 

be a single-plane time delay function induced by a lens potential ip with g > singularities, all of which 
are infinite singularities. Suppose that y is not on a caustic and Ty is isolated. Then for a generi^ Ty 
the number of lensed images obeys: 

1. N = 2N++g-l = 2iVsad -g + l, N+ = Ns^d -g + l- 

2. N>g+1, TV^in > 1, iVsad > iVmax + 9- 

3. If the corresponding lensing map is locally stabl^ with Crit(?7) bounded, then for sufficiently large |y| 
the lower bounds on the number of images are attained: 

N ^ g + l, TVinin = 1, A^max = 0, and A^sad = g- 

As long as the hypotheses of Theorem [T] hold, the image counting information is independent of the 
choice of gravitational lens model. The topological nature of the counting formulas and lower bounds 
give them wide applicability. Part 1 of Theorem[l]states that the number of images has parity (even-ness 
or odd-ness) opposite to the number of singularities, part 2 gives lower bounds — e.g., there are at least 
g+l images and at least g saddle images, and part 3 implies that the lower bounds in Part 2 are actually 
the smallest number of images that are achievable. 

Application to microlensing: Theorem [T] applies to subcritical microlensing, i.e., lensing due to g point 
masses with continuous matter Kc and shear 7 satisfying 1 — Kc + 7 > and 1 — Kc -f 7 > . There is 
an even number of images if and only if the number of stars is odd. Also, since there is no maximum 
images, we have 

Nmin = N,^d - g + 1- 

^ Theorem [1] holds either for Ty or almost all sufficiently small linear perturbations of Ty, i.e., the functions Ty (x) + p ■ x 
for every p g M^, except in a set of measure zero, with |p| sufficiently small |64l p. 421]. 
^ Local stability is equivalent to the set of critical points of tj consisting only of folds and cusps 1641 p. 294]. 
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This is a useful formula for checking whether images are overlooked by in microlensing simulations. 
Consult [Ml Chap. 1 1] for a detailed discussion, where the cases 1 — Kc + 7 > and 1 — Kc + 7 < (strong 
shear lensing) and 1 — Kc + 7 < and 1 — Kc + 7 < (supercritical lensing) are also treated. 

In the situation of a nonsingular lens, the following corollary of Theorem [Ijl) recovers (by setting 
g — 0) the single-plane Odd Number Image Theorem found in 1981 by Burke 16 , who proved the result 
using a different approach, namely, the Poincare-Hopf index theorem. 

Corollary 1 16J (Single-Plane Odd Number Image Theorem) For a non-caustic point y, let Ty 
be a single-plane isolated time delay function induced by a nonsingular gravitational lens potential. Then 
the total number of images is odd: 

N^2N+-1^ 2iVsad + 1. 

Nonsingular lensing is typically due to modeling galaxies as smooth on a macro scale. Though the 
predicted number of images is odd, typically an even number of images is observed. The reason is that 
maximum images are angularly located where the surface mass density of the lens is supercritical (core 
of galaxies), which causes them to become significantly demagnified (e.g., [51j, |64i P- 470]). 

We saw from Theorem[Tl2) that if an isolated lens has at least one singularity {g > 1), which can be 
a point mass, singular isothermal sphere, etc., then the lens can produce multiple images. Namely, there 
is a point in the light source plane from which a light source has more than one lensed image, TV > 1. 
How about multiple imaging due to an isolated nonsingular lens? The following necessary and sufficient 
condition for multiple imaging by such lenses was established in 1986 by Subramanian and Cowling [78] : 

Theorem 2 [75] (Criterion for Multiple Images: Single-Plane Nonsingular Case) For a non- 
caustic point y, let Ty be a single-plane isolated time delay function induced by a nonsingular gravitational 
lens potential. Then: 

1. N > 3 if and only if there is a point Xq in the lens plane L such that det[HessTy](xo) < 0. 

2. If there is a point Xq G L where the surface mass density is .supercritical, k(xo) > 1, then a light 
source at yg = ''7(xo) will have multiple images, N > 3. 

We prove Theorem [5] to illustrate how some of the previous counting results and ideas are used theoret- 
ically. 

Proof (1) If > 1, then Corollary (TJl) yields N — 2Nsa,d + 1 > 3, so iVgad > 1. Because there is at least 
one saddle image, say, Xgad, the point Xgad G L satisfies det[IIessry](xsad) < 0. Conversely, if there is a 
point Xq E L such that det[HessTy](xo) < 0, then a light source at yo = ■'7(xo) has at least one saddle 
image. Since A^sad > 1, -^min > 1, and the number of images is odd, we have iV > 3. 
(2) If k(xo) > 1, then image xq cannot be a minimum since minima are located where the surface mass 
density is subcritical (e.g., [64| p. 423]). Then xq is either a saddle or maximum. For Xq a maximum. 
Theorem [Tj2) with g = (nonsingular case) yields A^sad > .^max > 1, which yields TV > 3 (since 
^min > !)■ For Xq a Saddle, we get > 1 because there is also at least one minimum and the number 
of images must be odd. □ 

The sufficient condition k(xo) > 1 for multiple images of a source at yo — f]{:Ko) is not also a necessary 
condition. A perturbed Plummer lens has the following respective potential and surface mass density: 

V;(x) = ^ log(l + IxH - |(u2 - v^), «(x) = (Y^)^, 

where 7 > is the external shear and x = (u, v). This lens is subcritical everywhere for < kq < 1: but 
can still produce multiple images (e.g., [HH p. 429]). 

Multiple imaging is also discussed in Section [5] from the optical geometry point of view, and for 
nonsingular lenses in a spacetime context in Section [T) 
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2.2 Maximum Number of Images: Single Plane Lensing 



A natural next question is to determine the maximum attainable number of images due to g stars. Using 
a trick with complex quantities, Witt [52] showed in 1990 that g point masses will generate at most 
+ 1 images. Since it was unclear how to extend the trick in (55] to multiplane lensing, AOP gave 
in 1997 an alternative proof of the upper bound using resultants, an approach gcneralizable to multiple 
lens planes (see Section 15^ . 

Combining the above upper bound result with Theorem [Il2), the number of images is bounded as 
follows: 

g+l<N<g^ + \. (2) 

For two point masses g = 2, equation ^ yields = 3,4 or 5. However, we cannot have A' = 4 since 
by Theorem[T]Jl) an even number of stars has to produce an odd number of images. Hence, A = 3 or 5, 
which recovers the result in [73j . 

It remained unclear whether the quadratic upper bound oi g^ + 1 in ^ can be attained. In 1997, Mao, 
AOP, and Witt [46j conjectured that the maximum number of lensed images should be linear in g. They 
also constructed a lens system consisting of g point masses of equal mass I/5 on the vertices of a regular 
polygon and showed that this symmetrical system produces a maximum number of 3g + 1 images. The 
latter is a linear lower bound on the maximum number of images for g point masses. Rhie |70] showed 
in 2003 that by putting a mass > at the center of the regular polygon, but with equal masses 
l/((7 — 1) on its vertices, a total of 5g — 5 images can be achieved for sufficiently small masses. Using a 
modification of the method in [46] , Bayer and Dyer [13] gave in 2007 a much simpler proof of the result 
in [70] and improved our understanding of the result by determining an upper bound of the central 
mass such that the maximum number 5g — 5 of images is attained for all rrie < m^. The conjecture was 
finally settled in 2006 when Khavinson and Neumann [42j employed complex rational harmonic functions 
to show that the total number of images is at most 5g — 5. The results are summarized below: 

Theorem 3 (MELximum Number of Images: Single-Plane Case) For g > 2 point masses on a 
lens plane and a light source not on a caustic, the number of images satisfies: 

1. |42] Af < 5g- 5. 

2. The upper hound 5g ~ 5 is attainable for: 

(a) [46 g = 3 if the point masses have equal mass l/g and lie on the vertices of an equilateral triangle 
centered at the origin and inscribed in a circle of radius r bounded above as follows: 



r <rr 



(b) |701ll3j g > 4: if g — 1 of the point masses have equal mass l/{g — 1) and lie on the vertices of a 
regular polygon centered at the origin and inscribed in a circle of radius r < Tcr, and a point mass 
of rUe is at the center of the polygon with mass upper bound given by: 



< me < m^, 



2(1+1.) 



1 



We now have the following sharp (i.e., attainable) bounds on the total nmnber of images due to point 
masses: 

g+l<N <5g-5, g > 2. (3) 

Remark: For lensing by a general matter distribution, there is no overarching maximum number of 
images because a mass clump can always be added to a lens to create more images. 



3 Multiple Imaging in Multiplane Lensing 

To set up the image counting results for multiple lens planes, we review some of the notation and concepts 
needed. Consult jM] Chap. 6] for more details. 

Let Li be the ith lens plane counting from observer to the light source plane and set X — (xi, . . . , x^) S 
ii X • • • X Lk', see Figure[2] Denote the gravitational lens potential on Li by ipi, where j = 1, . . . , fc, and 
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Fig. 2 A schematic of fc-plane gravitational lensing. The action of the lensing map is captured by tracing Hght rays 
backwards from a subset P of Li to the light source plane S. Credits: After |64l p. 199]. 



the light source plane by S' = with its elements y G S. Let Ty : Li x • • • x L^. ^ R, be the fc-plane 
time delay function induced by the potentials ipi ^-^d denote the associated fc-planc lensing map by 
■ p ^ s Here P = R"^ — B C Li, where B is the set of light path obstruction points. Note that for 
single-plane lensing P = Li = L and B ~ A (set of singularities of the potential -01 = "0)- 

The set Crit(?7^'^^) of critical points of jy^*^-' is the locus of all infinitely magnified lensed images for all 
light source positions in S, while the set of caustics of jy^*"'' is Caustic(?7'^'^^) = rj'*^' [Crit(r7*^''))] C S, i.e., 
the set of all light source positions from which there is at least one infinitely magnified lensed image. 
When Crit(r7('^^) is bounded, then both Crit(r7('')) and Caustic(r7('^') are compact p. 293]. 

(k) 

The fc-plane time delay function Ty is called subcritical at infinity if for each non-caustic point y and 
|X| sufficiently large, the eigenvalues of Hess Ty'^'' (X) are positive, i.e., the time delay surface has positive 
Gauss-Kronecker curvature for |X| sufficiently large. If Ty*^'' is subcritical at infinity and ry'^''(X) oo 
as |X| oo, then Ty*^-* is called isolated. 

(k) 

By Fermat's principle, the critical points of Ty determine the lensed images of a source at y 
(e.g., [64], Chaps. 3,6), so the lensed images are then identified as generahzed saddles of index i, where 
i = 0,l,...,2fc. 

Notation: 

1. Ni = number of images of index i. Here A'o and N2k are, resp., the number of minima and maxima. 

2. = J2i{cvcn) ~ number of even index images. 

3. ~ X)i(odd) ~ number of odd index images. 

4. = A^_l_ + N- — total number of images. 



3.1 Counting Formulas and Lower Bounds in Multiplane Lensing 

In 1991, AOP [61J extended to fc-lens planes the image counting results of Theorem[Tl including the Odd 
Number Image Theorem. All the results are topological in nature, which give them wide applicability. 
As before, we state the results only for the subcritical case (see [64 ' Chap. 12] for more): 

Theorem 4 [HT] (Multiplane Counting Formulas and Lower Bounds) Let Ty'''' : ii x • • • x Lj. ^ M 

be a k-plane time delay function induced by lens potentials ij^i, . . . ,tpk, where each ipi has gi singularities, 
each of which is an infinite singularity. Assume that y is a non-caustic point and Ty^^ is isolated. Then 
for a generi^ Ty'^^ the number of lensed images satisfies: 

1. 2N+ - J](l - g,) = 2N^ + [](1 - g,), N+ = + [](1 - g,). 

k 

2. ^>n(i+.9*)- 

1=1 

3. No > I, Nj > 5£i • ■ • gij for I < j < k, and Nj > for k + I < j < 2k - I. 

l<ii<---<lj<k 



* Theorem holds either for Ty*"' or Ty'^^X) + p ■ X, for sufficiently small p G K^*', except in a set of measure zero 1641 
pp. 452-455]. 
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4- If the corresponding lensing map t]^^^ is locally stable and Crit(r7('^^) is bounded, then for sufficiently 
large |y| the above lower bounds on the number of images of different types are attained: 

k 

N = X{(l^g^), N^^\, Nj^ J2 9ir--9i, for l<J<k, 

1=1 i<ei<---<ij<k 
and Nj ^0 for fc + 1 < j < 2fc - 1. 

Corollary 2 (Multiplane Odd Number Image Theorem) For a non-caustic point y, let Ty^^ be a 
isolated k-plane time delay function induced by nonsingular gravitational lens potentials. Then: 

N = 2N+ - 1 = 27V_ + 1. 



3.2 An Upper Bound on the Number of Images 

A natural question is whether there are multiplane lensing bounds on the number of images analogous 
to the single-plane ones in ([2]) or Q. The following 1997 theorem of AOP [63] is what is known so far: 

Theorem 5 |63i (Multiplane Upper Bound) Let y be a non-caustic point. Then the total number 
of images due to gravitational lensing by g point masses on g lens planes with one point mass on each 
lens plane is bounded as follows: 

2^ <N <2 (2^(s-i) „ 1^ ^ g>2. (4) 

By Theorem 13^1,4), the lower bound of 2^ is sharp and the total number of images is always even. 
Theorem [5] was proven using the theory of resultants. 



4 General Stochastic Lensing: The Expected Number of Minimum Images 

Stochastic lensing occurs when a component of a systems is random, typically, the lens. This relates 
to the broader study of random functions, a subject that has been explored in mathematics primarily 
for Gaussian random fields (e.g., Adler and Taylor [5], Azais and Wschebor [T^], Forrester and Honner 
pS] . Li and Wei [44^, Shub and Smale [75], Sodin and Tsirelson [77], and references therein). However, in 
gravitational lensing most of the realistic lensing scenarios produce non-Gaussian random fields that have 
not been previously considered (e.g., see Theorem [7| below). Therefore, a new mathematical framework 
needs to be developed for the study of stochastic lensing. 

A natural first step in stochastic lensing is to study the expectation of the random number of lensed 
images. We shall present some recent rigorous mathematical work in that direction. 

The expectation of N+{D,y), the number of positive parity lensed images inside a closed disk D of 
a light source at position y, is given by the following Kac-Rice type formula (see 

E[N+{D,y)]= f i?[((l-«;(x))2-r2(x)) le,(x) r7(x) = yl /^(,) (y) dx, (5) 
Jd l J 

for almost all y. Here, Ig^ is the indicator function on C/^ = {i^ e : G(i^) € (0,oo)}, where G'(x) = 
det[JacJ7](x), and /r,(x) is the probability density function (p.d.f.) of the lensing map at x. 

This formula holds for a fixed light source position. Unfortunately, this position is unknown in most 
lensing observations. Therefore, a physically relevant extension of ^ is to view the light source position 
as a random variable whose probability density function is compactly supported over a subset of the 
light source plane having positive measure. This result can be further generalized to the entire light 
source plane. To do so, a countable compact covering {S} of light source plane is considered, with each 
set 6 in the family having the same positive area. A corresponding family {Ye} of random light source 
positions Y is constructed, with each Y uniformly distributed over a set 6. The global expectation of the 
number of positive parity lensed images, denoted by E[N+{D,Y; S)]{g]., is then defined as the average 
of E[N+{D,y)] over the family {Ys} (equivalently, over the family &). This notion was introduced in 

m- 

AOP, Rider, and Teguia 2009 (BBj determined a general formula for the global expected number of 
positive parity images: 
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Theorem 6 [53] (Global Expected Number of Positive Parity Images) The global expectation 
of the number of positive parity lensed images in D is given by: 

E[N+{D,Y-6)]{e}^T^^ / [det [Jacr?] (x) le,(x)] dx, (6) 

where \&o\ = area(S). 

The importance of Theorem|6]is its applicability to a wide range of lensing scenarios, with few assumptions 
on the distribution of the random gravitational field within each scenario. Below we will discuss an 
application of this theorem to image counting in microlensing. 



5 Stochastic Microlensing: Asymptotics 

The mathematical analysis of the p.d.f.s and expectations in microlensing is very difficult. An asymptotic 
approach is a natural first step that will be important for understanding first order terms and how 
deviations from them occur at subsequent orders. For example, we shall see that, though the lensing 
map of microlensing is a bivariate Gaussian at first order, the mapping deviates from Gaussianity at 
the next orders. The p.d.f.s of the time delay function and shear will also be seen to be non-Gaussian. 
Unfortunately, this means that most of the technology already developed in the mathematical theory of 
random fields will not be applicable directly to stochastic lensing. 



5.1 Notation 

The stochastic microlensing scenario we shall discuss is one with uniformly distributed random star 
positions. Recall from the Introduction that the potential is given by 

^^(x) = ^|xp - ^(u^ - v^) + ^m, log |x - 
where x = (u, v) , the time delay function at y by 

g 

Tg,y (x) = di (x; y) - ^ mj log |x - \ , 

where c?i(x;y) = i|x — yp — -^1^ |xp + ^(u^ — v^), and the components of the lensing map 77^ — {'nig,T]2g) 
by: 

where = {Uj,Vj). 
Notation and Assumptions: 

1. Equal masses: rrij = m, where j — 1, ■ ■ ■ g. 

2. R = \/g/TT and — Tim. 

3. S(0, R): closed disc of radius R centered at the origin 0. 

4. The random point mass positions ^1, . . . , ^re independent and uniformly distributed over -6(0, R). 



(1-Kc + 7)U-X! 



[/,-u 



^ (t/,-u)2 + (y,-v)2' 



(1-Kc-7)V-Xl 
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Now, normalize the random time delay function and random lensing map as follows: 
T;y(x) ^ T3^y(x) + gm logi?, r,^^) ^ 



Write the components of the random shear tensor due only to stars as follows: 

r. ,.S- rn[{U,~uY-{V,-Yf\ . . 2m{U , - u){V, - v) 

- [u, - u)2 + (y, - v)2]2 ' ^'^^"^ - [u, - u)2 + [v, - vf] 



For fixed y and x, we denote the possible values of the previous random quantities as follows: 

1. t = possible values of r^y(x). 

2. (h, k) = possible values of T7g(x). 

3. (z,w) = possible values of (A,g(x), r2,g(x)). 

4. G = (z^ + w^)^/^ = possible values of the magnitude of the shear. 

For the asymptotic p.d.f. of the normalized random lensing map, we shall need the quantities 



(l-«:, + 7)u (l-K,-7)v , /log(B5i/2) 2^e^-^^ 
ai = 1== , 3.2 = == , a , B , 

where 7e is Euler's constant. For the asymptotic p.d.f. of the random shear, we shall use: 

A.,|xp kUQkI - 9G^) KljSKt - 2Anl(f + 3G^) 15«t {^4 - 40^^G^ + ISG^) 
' '""'^ m 2(k2 + g2)2 4(^2+g2)3 + 32(k2+g2)4 • 

Note that H2 depends on |x|2/rn, which is the mean number of point masses within the disc of radius 
|x| centered at the origin. 



and 



5.2 Asymptotic P.D.F. s of Random Time Delay Function, Random Lensing Map, and Random Shear 



In 2009, AOP, Rider, and Teguia j65P66j used a rigorous mathematical approach to characterize up to 
order three the asymptotic p.d.f.s of the microlensing normalized random time delay function, normalized 
random lensing map, and random shear: 

Theorem 7 Fix y e S and x = (u, v) e B{0, R). Then: 

1. [SS] (Random Normalized Time Delay Function) In the large g limit, the asymptotic p.d.f. of 
T* y{x) takes the following form: 



(x)(t) 



(t-di-e)«-i 
(9-1)! 



exp 



2(t-di-c) 



0, 



t > di + c 
t < di + c 



+ 0{g 



-3/2\ 



The first term of the p. d.f. is a Gamma distribution. 
2. [S^ (Random Normalized Lensing Map) In the large g limit, the asymptotic p.d.f. of t]*{x) 

takes the following form: 



(li-ai)^+(k-a2)^ 



f77*(x)(h,k) = 



2TTa: 



1 — K: 



(h — ai)u + (k — a2)v 



^ ((h~ai)2 + (k-a2)^-2(ap2 

'4^ (KT 



log(log5) 



V Off (7/ 



The first term of the p.d.f. is a bivariate Gaussian distribution, but the next two terms highlight that 
J7*(x) already becomes non- Gaussian for large finite g. 
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3. [66] (Random Shear) In the large g limit, the asymptotic p.d.f. of (r']*g(x), /^2*g(^)) to-kes the 
following form: 



fA,.,(x).r2,g(x)(z,w) - 



27r(K2 + G2)3/2 



1 



Hi{G) , H2{G; \x\ 



9 



9' 



+ 0(g-3), (10) 



where G = (z^ + w^)^/^ denotes the possible values of the magnitude of the shear. The first term of 
the p.d.f. is a stretched bivariate Cauchy distribution. 



By equation (|10p . the asymptotic p.d.f. of the magnitude of the shear, namely, -r'g(x) = ^ /'^^(x) + g{'x.), 
is given as foUows (66] : 



fr,(x)(G) = 



(k2+G2)3/2 



Hi(G) , H2{G; \x\ 



9 



9^ 



0(9-')- 



Remarks: 

1. The first term in ([9]) was basically found in 1986 by Katz, Balbus, and Paczynski ^38], who actually 
determined the first term of the p.d.f. of the bending angle due only to stars, namely, the p.d.f. of 
the random vector ag{x) = T7g(x) — ((1 — Kc + 7)u, (1 — Kc — 7)v). 

2. The first term in (fTO|l was found in 1984 by Nityananda and Ostriker [52] . 

Subsequent work in 2009 by Keeton [39] used semi-analytical and numerical methods to study the 
stochastic properties of the lens potential, deflection angle, and shear under different assumptions about 
the distributions of the stars' masses and positions. 



5.3 Global Expected Number of Micro-Minima 

Wambsganss, Witt, and Schneider 1992 [79] determined the limit (/ ^ oo of the (globalfl) expected 
number Eq of minima in the entire plane R2 for microlensing without shear. This was extended to the 
case with shear in 2003 by Granot, Schechter, and Wambsganss [35] : 

(l-..)^-(7 + z)2-^^^^^^ 



" 27r|(l-.Jtot)2-72|76 (^2+G2)3/2 

where Ktot — + Kc and 7 is same shear employed at the microlensing scale. 

The value Eq resulting from the limit g — > 00 is, of course, independent of g and can be treated as the 
first term in an asymptotic expansion in 1/g. On the other hand, since the first term Eq is independent 
of 5, if we have only the term Eq, then there is no way to know analytically the smallest g needed to 
lie within a certain approximation of the global expected number of micro-minima. This is especially 
important in numerical simulations. Therefore, we need to find more terms in the asymptotic expansion. 

AOP, Rider and Teguia 2009 [BH] applied Theorem [H] and Theorem [71^3) to determine rigorously the 
global expected number of micro-minima up to three orders: 

Theorem 8 [SB] (Global Expected Number of Micro-Minima) Let D be a closed disc and suppose 
that continuous matter is subcritical, i.e., < < 1- Then, in the large g limit, the first three asymptotic 
terms of the global expectation of the number of minimum images in D is given by}^ 

^\AT m V «M ^* Mg,6o f (l-Ate)"-(7 + z)2-w2 [ HijG) , g2(G;ao) ]^ ^ , „3, 
E[Ng^^,r,{D, Y; 6)]{6} = ^ (^2 ^ ^2)3/2 + + ^ J ^^^w + 0{g ), 

(11) 

where fJ-D,eo — J§}\' ^ is a closed disc of radius 1 — Kc centered at (—7, 0), and oq ~ Jjj ju^du. 

^ The terminology "global expected number of minima" was actually introduced in 1661 . where the difference between 
the expected number of minima and global expected number of minima was clarified. 
^ The quantities Hi and H2 were defined in ^ and JSj. 
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Theorem[5]gives us the leading three terms in an asymptotic expansion of the global expected number 
of micro-minima in any reference disc Z), not necessarily in the entire plane. The first term in that 
expansion (jlip is more general than Eq because Eg applies to the whole plane. A possible physical model 
for the factor ^D,eo in (HH) is the macro-scale magnification given by: 

1 

|(l-Ktot)^-7^l 

It is only in this case that the first term in (fTTjl coincides with Eq. 

Another consequence of Theorem [8] is it enables us to estimate how small we can choose g in order to 
have the first three terms in the asymptotic expansion (jlip lie within a certain percentage of the global 
expected number of micro- minima. Theorem [8] also shows us analytically how accurate an approximation 
the first term is to the exact global expectation by quantifying how the next two higher-order terms 
perturb the first one. 

Remark: The third-order term in depends on the mean number of point masses within the disc of 
radius |x| centered at the origin, namely, the quantity k,|xP/to; compare equation ([5]). 



6 Multiple Images in Optical Geometry 

6.1 The Optical Metric and Fermat's Principle 

The trajectories of spatial light rays can also be studied in optical geometry, which is conceptually 
between the thin- lens, weak-defiection approximation used in the previous sections and the full spacetime 
treatment of null geodesies. Optical geometry, which is also known as Fermat geometry or optical reference 
geometry, is a useful tool to investigate gravitational and inertial forces in General Relativity ^ . Recently, 
it has also been applied to field theory near black hole horizons by Gibbons and Warnick [3^, and to 
analogue models of gravity in Finsler geometry by Gibbons, Herdeiro, Warnick and MCW 

Fermat's principle and optical geometry in the conformally stationary case was discussed in detail by 
Perlick [56j . To illustrate this approach, consider for simplicity a static spacetime metric 

ds^ = ffoo (dt)^ + gij dx'dx^ , 

where the summation convention is used, and a null curve parametrized by < < 1, say, with tangent 
vector k. Then the variational principle yields (e.g., Frankel [23]): 

1 

g(^,k) (1) Jo 

if the null curve is in fact a null geodesic with Vkk = 0. Hence, one obtains Fermat's principle of 
stationary arrival time, rather than stationary travel time as in the case of flat space used in the impulse 
approximation. Now recasting the spacetime line element as follows, we see that the spatial light rays 
are geodesies of the optical metric with Riemannian signature, 

{dtf - dx'dx' = g°f' dx'dx', 

900 

by Fermat's Principle. 



6.2 Gauss-Bonnet and Lensed Images 

Optical geometry offers a different perspective on image multiplicity in gravitational lensing, which is 
also partially topological. This can be understood with the Gauss-Bonnet Theorem, which connects the 
local optical geometry with global properties of the light rays. Gibbons [IS] discussed this for cosmic 
strings, and this method was recently extended to spherically symmetric metrics by Gibbons and MCW 
[5T] and Gibbons and Warnick [3D] . In this approach, the occurrence of two images can be expressed as 
a digon of light rays: 
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Fig. 3 Isometric embedding in K'^ of the equatorial plane of the Schwarzschild solution in the optical metric. The waist 
occurs at the radius of the photon sphere, r = 3m. The Gaussian curvature is negative everywhere. 



Theorem 9 j281l31l[5U] (Light Ray Digons) Let {S, 5°^*) be a totally geodesic, simply connected surface 
endowed with an optical metric (7°^' and Gaussian curvature K . Let D <Z S be a geodesic digo"^ hounded 
by two light rays intersecting at the light source s S and the observer £ S with corresponding positive 
interior angles 9s and Oo- Then: 

Example: The Schwarzschild solution with mass parameter m is an instructive example. The equatorial 
plane in the optical metric outside the photon sphere is shown in Fig. [31 and it can be seen that the 
Gaussian curvature is negative everywhere. Indeed, a calculation of the Gaussian curvature gives 

^I^A-^UO. 



r3 (1 _ 2mj3/2 y 2r 

Spatial light rays, which are the geodesies of S, must therefore diverge locally, and the equation of the 
Theorem cannot be fulfilled. So the fact that two light rays can intersect at the source s and the observer 
o in Schwarzschild geometry shows that S cannot be simply connected, which indeed it is not because 
of the event horizon at r = 2m. Hence, the topological contribution to the Gauss-Bonnet theorem turns 
out to be essential for image multiplicity. 

Example: Now consider the Plummer model with mass toq a-nd scale radius tq. This is a reasonable 
model for an extended non-relativistic gravitational lens (i.e., a galaxy) which allows multiple imaging: 
three images are produced if the light source is in the maximal caustic domain. Of course, there is no 
event horizon in this model, so the surface S is simply connected, and the optical metric approaches that 
of the Schwarzschild solution with < at large radii r, since mo is finite. Therefore, the equation of 
the theorem holds, and we expect that the Gaussian curvature of the optical metric changes sign in order 
to allow for multiple images. A calculation of K to lowest order yields 



2mo 



^0 ^ ^ \ l + (r/ro) J 

which confirms that -R" > for small radii. 

A more detailed analysis shows that the Gauss-Bonnet theorem can also be used to calculate deflection 
angles. In the case of the singular isothermal sphere, for instance, the optical metric describes a cone so 
that the gravitational deflection emerges as a consequence of the cone's deficit angle, rather like a spatial 
analogue of lensing by cosmic strings 31j. 



7 Multiple Images in Spacetime 

7.1 Necessary and Sufficient Conditions for Multiple Images 

The theories for spatial light rays provide useful frameworks for gravitational lensing, especially since the 
impulse approximation can easily be applied to models of astrophysical interest. However, the fundamen- 

^ A geodesic digon in a lensing scenario is a polygon with two vertices, bounded by the two geodesies intersecting at the 
source and the observer. 
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tal arena of optics in General Relativity is of course spacetime. The study of wavefront singularities in 
spacetime was pioneered by Friedrich and Stewart [24l, who classified stable singularities in Minkowski 
space and considered the relationship with the initial value problem in General Relativity. The general 
form of Fermat's Principle in spacetime was investigated by Kovner , and a precise proof for arbitrary 
Lorentzian manifolds {M,g) was first obtained by Perlick in 1990. This also laid the foundation for 
a rigorous study of the conditions for multiple images in spacetime, which made an earlier result by Pad- 
manabhan and Subramanian [54] more precise. Perlick showed that the existence of a conjugate point or 
a cut point along a null geodesic is a sufficient condition, which becomes necessary if further conditions 
on the topology or causal structure of M are imposed: 

Theorem 10 [SS] (Multiple Images in Spacetime) Let (M, g) be a four- dimensional time- oriented 
Lorentzian manifold. 

1. Sufficient Conditions: Let X be a future-pointing null geodesic affinely parametrized by s, and fix 
A(si) = pe M. 

(a) If there is a S2 > si such that A(s2) is conjugate to A(si) along X, then there is a timelike curve 
through each A(s), s > S2, which can be reached from p along another future-pointing null geodesic. 

(b) If there is a S2 > si such that A(s2) is the future cut point of p along A, then there is a timelike 
curve through each A(s),s > S2, which can be reached from p along another future-pointing null 
geodesic. 

2. Necessary Conditions: Fix a timelike curve 7 and a point p £ M. 

(a) If there are two future-pointing null geodesies from p to ^ which are null homotopic, then there is 
a future-pointing null geodesic from p to j which contains a point conjugate to p. 

(b) If {M,g) is strongly causal and if there are two future-pointing null geodesies from p to 7, then 
the intersection with 7 comes on or after the future cut point of p along at least one of the null 
geodesies. 

Having established conditions for existence, one can now proceed with counting results for null geodesies. 
7.2 The Odd Number Theorem in a Spacetime Setting 

The proofs of the Odd Number Theorem for single and multiple lens planes discussed in Sections [5] and [3] 
employ finite dimensional Morse theory in the impulse approximation. Interestingly, this result can also 
be extended to spacetime under certain assumptions. McKenzie 1985 showed this using the degree of a 
map between two spatial spheres, which requires stationarity of the spacetime, and gave a proof using 
Uhlenbeck's Morse theory for null geodesies on globally hyperbolic Lorentzian manifolds 08]. However, 
these conditions appear to be too restrictive for realistic spacetimes, as discussed by Gottlieb 1994 |34j . 
An extension to infinite dimensional Morse theory on the Hilbert space of null curves was developed 
Perlick 1995 ^7] and Giannoni, Masiello and Piccione 1998 [27]. The resulting Morse relations have been 
applied by Giannoni and Lombardi [26j to prove a version of the Odd Number Theorem. In 2001, Perlick 
[59] defined the concept of a simple lensing neighborhood, which formalizes a physically meaningful lensing 
geometry to avoid some of the technical complexities, and proved the following spacetime version of the 
Odd Number Theorem: 

Theorem 11 |59] (Odd Number of Null Geodesies) [59^ Let U be a simple lensing neighborhood in 
a four- dimensional time-oriented Lorentzian manifold {M,g). Fix a point p £ U and a timelike curve 7 
in U such that it has no endpoints on the boundary dU . If j intersects neither p nor the caustic of the 
past light cone of p, then the number of past-pointing null geodesies from p to "f completely within U is 
finite and odd. 

A rigorous treatment of the optics in a spacetime setting can be found in Perlick [60] . 
8 Global Magnification Relations for Special Lens Models 

In gravitational lensing, when a source gives rise to multiple images, the magnifications of these images 
often obey certain relations. The simplest example of such a relation is provided by a single point-mass 



14 



lens. In this setting a source gives rise to two lensed images, and it can be shown that their signed 
magnifications always sum to unity: 

Ml + M2 = 1 (one point mass) 

(e.g., [M], P- 191). The surprising fact about this result is that it holds irrespectively of the mass of the 
point mass and its position on the lens plane. It is also a "global" relation, by which is meant that it 
involves all of the images of a given source and not merely a subset of them. Witt and Mao 1995 [53] 
generalized this result to a two point-mass lens. They showed that for a source lying anywhere inside the 
caustic curve, a region which gives rise to five lensed images (the maximum number in this case), the 
sum of the signed magnifications of these images is also unity: 

Ml + M2 + M3 + M4 + M5 = 1 (two point masses) . 

Like its predecessor, this relation is "global" and holds irrespectively of the lens's configuration, provided 
the source lies inside the region giving rise to the maximum number of images. Rhie 1997 [69j subsequently 
extended this result to N point masses. Dalai 1998 [TO] and Witt and Mao 2000 [51] then showed that 
other common lens models, such as singular isothermal spheres and ellipses (SISs and SIEs) and elliptical 
power-law potentials, possess similar global magnification relations, even when these models include 
shear. 

Dalai and Rabin 2001 [20] provided a residue approach that systematized and expanded on the above 
work. They began with the Euler trace formula, which they proved using residue calculus. The Euler 
trace formula identifies sums of magnifications as coefficients of certain coset polynomials; we shall return 
to this in Section [10. 31 below. The method in [20 i uses meromorphic differential forms in several complex 
variables. By considering a meromorphic 2-form on consisting of polynomials whose common zeros 
are the image positions of the lousing map, the authors were able to use the Global Residue Theorem, 
which states that on compact manifolds the sum of all the residues of a meromorphic form vanishes. This 
allows one to replace the common zeros of a form in C^, viewed as a subset of the compact manifold 
CP^, by minus the sum of residues at infinity in CP^. In this way magnification sums were transformed 
to a condition about the behavior of the lens equation at infinity. Their results are summarized in the 
following thereom: 

Theorem 12 20J The models listed below possess the following zeroth and first magnification moment 
relations: 



point masses 




1 


Zs Zj 


point masses + shear 




1 


Zs + IZs 


1 


_^2 


(1-7T 


SIE 




2 


2zs 


SIE + elliptical potential 




1 




SIS + shear 




2 


2{Zs + -IZs) 


1 


_^2 


(1-^2)2 


SIE + shear 




2 


2{Zs + 76^*""' Zs) 


1 


_^2 


(l_y)2 
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Notation: ruj is the mass of jth point mass, 7 is the shear with orientation and Zj and Zg are, 
respectively, the position of the jth lensed image and the position of the source using complex variables. 

Residue calculus methods were also used by Hunter and Evans 2001 [37], wherein magnifications of 
images were realized as residues of complex integrands. By Cauchy's theorem, sums of magnifications are 
then equivalent to a contour integral. This method was used to derive magnification relations for elliptical 
power-law potentials which expanded upon the work of 84J . As first shown in 84J , for an elliptical power 
law potential oc (u^ + v^g"^)''/^, where q is the ratio of the minor to major axes, the total signed 
magnification denoted by B is exactly B = 2/(2 — b) for the cases & = 0, 1; for other values, it is an 
approximation. The contour integral method used in [371l21j covered not only all cases when B is an 
integer, including second and third magnification moments, reciprocal moments, and random shear, but 
also cases with non-integer values of B as well. In 2002, Evans and Hunter [21] extended their results in 
[37j to include elliptical power-law potentials with a core radius, and calculated magnification invariants 
for subsets of the images with even and odd parities. 



9 Local Magnification Relations: Caustics up to Codimension 3 

9.1 Quantitative Fold and Cusp Magnification Relations: Single Plane Lensing 

All of the above magnification relations are "global" because they involve all the lensed images of a given 
source. But they are not universal because their relations were derived in the context of specific types 
of lens models (point-mass lenses, SIEs, etc.) There is another type of magnification relation, a so-called 
"local" magnification relation, that is universal in the sense that it holds for a generic family of lens 
models. It is called a "local" relation because it holds for a subset of the total number of lensed images. 
Such relations arise when the source lies close to a caustic singularity. The two simplest types of such 
singularities are the fold and the cusp. For a source near a fold, there will be two images straddling the 
critical curve, while for a source near a cusp, there will be a triplet of images. Interestingly, the signed 
magnifications of this doublet and triplet always sum to zero (e.g., Blandford and Narayan 1986 [T5] . 
Schneider and Weiss 1992 [74], and Zakharov 1999 

A*i + /i2 = (fold), ^1 + /i2 -I- ^3 = (cusp). (12) 

These relations are important in gravitational lensing because they can be used to detect dark sub- 
structure in galaxies with "anomalous" flux ratios. These anomalies arise as follows. For quasars with 
four lensed images, it is often the case that the smooth mass densities used to model the galaxy lens 
reproduce the number and positions of the lensed images, but fail to reproduce the image flux ratios. 
Mao and Schneider 1998 [37] showed that in such situations the cusp magnification relation (fT^ fails. 
They attributed this failure to the assumption of smoothness in the galaxy lens, and argued that this 
smoothness breaks down on the scale of the image separation. This suggests the presence of substructure 
in the galaxy lens. The possibility of this became even more intriguing when Metcalf and Madau 2001 
[49] and Chiba 2002 showed that dark matter would be a natural candidate for this substructure. 

Being able to identify which anomalous lenses have substructure now became a priority. Keeton, 
Gaudi and AOP in 2003 [40] and 2005 [41] developed a rigorous framework by which to do so, using the 
observable quantities: 

„ ^ A^i + _ Fi - F2 ^ Ml + /i2 + M3 _ Fi - F2 + F3 

" iMil + IM2I " + F2 ' ^"^P " iMil + IM2I + ImsI" i^i + + F3 ' 

where Fi is the observable fiux of image i. The importance of these quantities is as follows. If a source lies 
sufficiently close to a fold or cusp caustic, then (fT2]l predicts that i?toid and i?cusp should vanish. If -Rfoid 
and i?cusp deviate significantly from zero (Monte Carlo methods were used to determine what constituted 
a significant deviation), then that would indicate the presence of substructure in that particular lens. 
On this basis it was shown in [IDIHT] that for the multiply imaged lens systems they analyzed, 5 of the 
12 fold-image lenses and 3 of the 4 cusp-image lenses showed evidence of substructure. 
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9.2 Quantitative Elliptic and Hyperbolic Umbilics' Magnification Relations: Single Plane Lcnsing 

Consider a family of time delay functions T^y which induces a corresponding family of lensing maps 
r]^. Here y is the source position on the source plane and the parameter c can be any physical input, 
such as the core radius or external shear0 Using rigid coordinate transformations and Taylor expansions, 
the universal, quantitative form of the lensing map can be derived in a neighborhood of a caustic (see 
[72ll64] ). As mentioned above, the quantitative forms of lensing maps near fold and cusp caustics obey 
the fold and cusp magnification relations (|12p . Aazami and AOP 2009 p] showed recently that the 
quantitative forms of lensing maps near certain higher-order caustic singularities, namely the elliptic and 
hyperbolic umbilics, also satisfy magnification relations analogous to (|12p. Their work is summarized by 
the following theorem: 

Theorem 13 1 For any of the smooth generic family of time delay functions Tc.y corresponding to the 
elliptic umbilic or hyperbolic umbilic caustic singularities, and for any source position in the indicated 
region, the following results hold: 

1. (Elliptic Umbilic) satisfies the following magnification relation in its four-image region: 

fil + /i2 + M3 + M4 = 0. 

2. D'l (Hyperbolic Umbilic) satisfies the following magnification relation in its four-image region: 

Ml + + M3 + M4 = 0. 

An application of this theorem to substructure studies was also given in [T|, using the hyperbolic 
umbilic {D'^) in particular. Analogous to the observables i?toid and ^cuspj the authors considered the 
following quantity: 

^ ^ A^i + + M3 + M4 _ f 1 - ^2 + f 3 - Fj 
|mi| + IM2| + IM3| + IM4| ^ Fi+F2 + F3 + Fi 

By Theorem I13[ -Rh.u. should vanish for a source lying sufficiently close to a hyperbolic umbilic caustic 
singularity and lying in the four-image region. One advantage of i?h.u. is that it incorporates a larger 
number of images than i?foid and i?cusp, and also applies to image configurations that cannot be classified 
as fold doublets or cusp triplets. In fact recent work has shown that higher-order caustics like the 
hyperbolic umbilic can be exhibited by lens galaxies. Evans and Witt 2001 [22 , Shin and Evans 2007 
[75], and Orban de Xivry and Marshall 2009 [53] have shown that realistic lens models can exhibit 
swallowtail (A5) and butterfly caustics (A4), as well as elliptic umbilics {D^) and hyperbolic umbilics 
(D^). It is hoped that such lensing effects will be seen by the Large Synoptic Survey Telescope, and thus 
that higher-order relations such as i?h.u. will become applicable in the near future. 

An example of the multiple imaging, critical curves, and caustic curves due to a hyperbolic umbilic 
(D^) is shown in Figured 

9.3 Universal Local Magnification Relations: Generic Caustics up to Codimension 3 
9.3.1 Preliminaries: Lensed Images and Magnification for Generic Mappings 

Consider a smooth, n-parameter family _Fc,s(x) of functions on an open subset of that induces a 
smooth [n — 2)-parameter family of mappings fc(x) between planes {n> 2). The functions Fc_s are used 
to construct a Lagrangian submanifold that is projected into the n-dimensional space {c, s} = M"~^ x M^; 
the projection itself is called a Lagrangian map. The critical values of this projection will then comprise 
the caustics of fc (e.g., Golubitsky and Guillemin 1973 [33], Majthay 1985 45 , Castrigiano and Hayes 
1993 [17 , and 64, pp. 276-86]). Arnold classified all stable simple Lagrangian map-germs of n-dimensional 
Lagrangian submanifolds by their generating family F^^s ([SJ, Arnold, Gusein-Zade, and Varchenko 1 1985 
[g p. 330-31], and [64, p. 282]). 

Given fc(x) = s, we call x e a lensed image of the source (or target) point s S R^; note that our 
lensed images are never complex- valued. Equivalently, lensed images are critical points of Fc^s, relative 

* There is at most one universal unfolding parameter c for caustics up to codimension 3. 
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Fig. 4 Multiple images of a source inside a cross section of a hyperbolic umbilic caustic in lensing. In each panel, the 
figure on the left depicts the caustic curve with source position (solid box) in the light source plane, while the figure on 
the right shows the critical curve with corresponding image positions (solid boxes) in the lens plane. The + sign indicates 
that only minimum and maximum images lie in the given region and the — sign labels a region where only saddle images 
occur. Credits: Figure from |T]. 



to a gradient in x. Next, we define the signed magnification 2Jl(xi; s) at a critical point of F^^s to be 
the reciprocal of the Gaussian curvature at the point (x.^, i^c.s(xi)) hi the graph of F^^s'- 

OT(x,;s) ^ 



Gauss(xi, _Fc,s(xi)) 



The advantage of this definition is it makes clear that the signed magnification invariants are geo- 
metric invariants. The signed magnification is expressible in terms of fc since Gauss(xi, i^c.s(xi)) — 
det[Hessi^c,s](xi) and each fc satisfies det[Jacfc] = det[HessFc_s]. It follows that 



det[Jacfc](xi 



which is the more common definition of magnification. If det[Jacfc](x) — 0, then x is called a critical 
point of fc- The collection of such points form curves called critical curves. The target fc(x) of a critical 
point x is called a caustic point. Though these typically form curves as well, they could also be isolated 
points. Since there are caustic curves for each value of the parameters c, varying these parameters traces 
out a caustic surface, called a hig caustic, in the larger space {s,c} — E". An example of critical and 
caustic curves is shown in Figure |4] for the hyperbolic umbilic (D^), along with various source and image 
configurations. 



9.3.2 Universal Magnification Relations for Generic Caustics Up to Codimension 3 

The following theorem about magnification relations for generic caustics up to codimension 3 was estab- 
lished by Aazami and AOP in 2009 [I]: 

Theorem 14 [I] For any of the smooth generic families of functions Fc.s {or the induced general map- 
pings fc) giving rise to a caustic of codimension up to 3, and for s, any non-caustic point in the n-image 
region, the following holds for the image magnifications dJli = 97l(xi ; s) ; 

n 

^971, = 0, 

1=1 

where n = 2 for a fold, n = 3 for a cusp, n = A for a swallowtail, elliptic umbilic, or hyperbolic umbilic. 
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The proofs of Theorems [T3] and [H] in [T] , including the subsequent shorter proof in [5] , are algebraic 
in nature and the method will be highlighted in Section [10.31 

After [T] appeared, an alternative proof was given by MCW 2009 [5T] that introduced new Lefschetz 
fixed point technology in gravitational lensing, clarifying an earlier approach 80J. The next section 
overviews the method in [81j . 



9.4 A Lefschetz Fixed Point Approach to Theorem [T51 and Theorem [HI 

We first review some needed basics from Lefschetz fixed point theory and then outline how it can be 
used to prove the aforementioned theorems. 

9.4-1 Holomorphic Lefschetz fixed point theory. 

If / : A/ ^ M is a smooth map on a compact manifold M, then its fixed points are fix(/) = {x G M : 
f{x) = x}, that is, the intersection of the graph {(x, /(cc))} S MxMwiththe diagonal {(cc, a;)} £ MxM. 
Fixed point theory, then, connects local properties of the fixed points, called fixed point indices, with 
global properties of / and M. In the case of a real manifold M, this is called the Lefschetz number L{f), 
which is a homotopy invariant because / induces a map on the space of closed forms and hence on the 
cohomology classes of M. For complex M and holomorphic /, the relationship between the analogous 
holomorphic Lefschetz number Lhoi(/) and the local fixed point indices is called holomorphic Lefschetz 
fixed point formula. The Lefschetz fixed point formulas are well-defined provided that the intersections 
are transversal, and can be regarded as special cases of the Atiyah-Bott Theorem [TOl[TT] . 

To illustrate this concept, we now discuss polynomial maps on the Riemann sphere C = CP^ = 
C U {oo} as an instructive example. Here, the holomorphic Lefschetz fixed point formula is also known 
as the Rational Fixed Point Theorem, which has important applications in complex dynamics (see, for 
example, the discussion by Milnor [5U]): 

Theorem 15 (Rational Fixed Point Theorem) Let / : C ^ C fee a rational map which is not the 
identity. Then: 

zefix(/) dz 

Here, the holomorphic Lefschetz number for complex projective space is ihoi(/) = 1. To see why this 
is true, recall that a rational map on the Riemann sphere can always be extended to a holomorphic 
map; then we can choose local holomorphic coordinates so that z = for some fixed point and write 
f{z) — §7(0)2; + O(z^) by holomorphicity. Hence, 

1 1 r dz 



(l-^)(0) 27ri7c 

where C is a loop enclosing only the fixed point at the origin. We may take / to be bounded so that 

z-f{z) z z{z-f{z)) z^ 

Using these results, if we now extend the loop C to Coo of infinite radius enclosing all fixed points of /, 
then 

1 1 / dz 1 / dz 



^ 2m fc z~f{z) 2m fc z ^ 

as required, since the integral of (fT3)l vanishes [50] . 

For a general complex manifold M with dimension d, the holomorphic Lefschetz formula is (see 
Griffiths and Harris ^36] , for instance) 

^-(^)= ^ d^^iiT^r ^''^ 

zefix(/) ^ ^' 

where is the d-dimensional identity matrix and Df is the matrix of first derivatives with respect to 
local holomorphic coordinates. Again, the transversality condition can be expressed as the requirement 
that the fixed point indices be well-defined, that is, det[Ld — Df] ^ 0. 
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Singularity 




deg{/C),dcg(/2C) 


n 


Fold 


21, 2| 


1 9 


2 


Cusp 


21, 2122 + 23 


1,3 


3 


Swallowtail 


2122 + C2f + 2J, 22 


4,1 


4 


Elliptic umbilic 


322 — 32 J — 2c2i, 62122 — 2C22 


2,2 


4 


Hyperbolic umbilic 


— 3zf — C22, — 32I — C21 


2,2 


4 


Elliptic umbilic (lensing map) 
Hyperbolic umbilic (lensing map) 


2? - 2|, -22:12:2 + 4C22 
2 J + 2C22, 2I + 2C21 


2,2 
2,2 


4 
4 



Table 1 Components of the complex maps fj(2i, 22) for generic singularities up to codimension three, their degrees and 
the numbers n of real solutions in the maximal caustic domains at finite positions. Here c denotes a control parameter. 
The lower section lists the corresponding properties for the quantitative elliptic and hyperbolic umbilics for a lensing map. 



9.4-2 Application to Theorems \13\ and \14\ 

We now outline the Lefschetz fixed point proof of Theorem [TH In the present approach (compare Atiyah 
and Bott |JJJ), all real solutions fc(a;i,a;2) = (si,S2) are treated as the real fixed points of a suitable 
complex map. So, first, we need to find a complexification that allows the application of the holomor- 
phic Lefschetz fixed point formula. The standard complexification {xi,X2) '—^Xi+ 1x2 does not yield 
holomorphic maps, but this problem can be circumvented, at the expense of the dimension, by treat- 
ing {xi,X2) = (zi,Z2) as independent complex variables on C'^. The corresponding complex generic 
maps = (/fj/l') : ^ are shown in Table [l] and the maximum number of solutions of 
/f (2:1,2:2) = si, /2'(zi, 22) — 52, possibly complex, is deg(/f ) deg(/2') by Bezout's Theorem. Now, since 
Table [1] shows that this is always equal to the maximum number n of real and finite solutions, we see 
that our complex formalism gives the usual real solutions in the maximal caustic domain, as required. 
Next, one can define the map / = (/i, /2) = (21 — /f + 2/1 , 22 — /2' + 1/2) on such that its fixed points 
are in fact those solutions. Hence, by construction, / is holomorphic and has no fixed points at infinity. 
Notice also that, at these fixed points. 



Fo 


— ^0 ' 




Fi 


V ryJTL — l 

- ^1^0 


^ ^0 


F2 







det[J2-i:'/] det[Jacfc] ' ^ ' 

Although this looks rather suggestive, we cannot apply the holomorphic Lefschetz fixed point formula di- 
rectly to /, since is not compact. However, one can rewrite / in homogeneous coordinates {Zq, Zi, Z2), 
where 21 = Zi/Zq and 22 — Z2/ZQ for Zq ^ 0, and consider the a map F — (Fq : Fi : ^"2) on CP^, which 
is of course compact, where 



z;p-'^^^(f?)fC(^Zo,Z,,Z2) + y,ZS\ 
'"^f^iZo,Z^,Z2)+y2Z^, 

and m = max(deg(/f ), deg(/2^)) > 2. Thinking of CP^ = U CP^, with : Zq = 1, CP^ : = 0, we 
recover F\c2 = f. Also, F is holomorphic since / has no fixed points at infinity. Finally, the holomorphic 
Lefschetz fixed point formula p4)) applies to F and to i^lcpi, where Lhoi = 1 as shown in Section fQ. 4. II 
Hence, using l|15[) . we obtain 1 — X^ILi + 1- This shows that the universal magnification invariant is 
zero for the generic cases in Theorem I14[ which can now be interpreted as a difference of two Lefschetz 
numbers. 

For the quantitative case a similar argument yields a proof of Thcorem[T3l the lower section of Table[l] 
shows that the properties of the quantitative elliptic and hyperbolic umbilics for lensing maps allow a 
construction as discussed in Section [9A2l in exactly the same way as for the corresponding generic maps. 

Remark: The above fixed point approach to lensing should not be confused with the one studied by 
AOP and Wicklin 1998 [ST], where fixed points of the lensing map were explored to determine source 
positions that have a lensed image coinciding with the unlensed angular position of the source. It is also 
interesting to note that the technical transversality condition mentioned in Section [9 . 4 . 1 1 simplv becomes 
< c» in the lensing context, that is, the usual condition for regular images. 
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Class Fc,s{xi,X2) 

An ±x"^^ ± a^i + Cn-lXi^^ + ■ ■ ■ + Cgxf + S2xj - SlXi ± S2X2 

(n > 2) (±{n + l)xY + {n - l)cr,-ix'^-^ + ■ ■ ■ + 3c3xj T 4:'-^2Xi , T 2x2) 

Un X^X2±X2 + Cn-2X2 + ■ ■ ■ + C2X2 — S2X2 — S\X\ 

(n > 4) (2x1x2 , xf ± (n - l)x^~2 + (n - 2)c„_23;2"^ + ■ ■ ■ + 2C2X2) 

Eq xf ± Xj + C3XIX2 + C2X| + C1X1X2 — S2X2 — sixi 

(sx^ + C3x| + C1X2 , ± 4x2 + 2C3X1X2 + 2c2X2 + cixi^ 

Ej X'l + X1X2 + C4X2 + C3X2 + C2x| + C1X1X2 — S23;'2 — •5lXi 

(sxi + Xj + C1X2 , 3xix| + 4C4X2 + 3c3x| + 2c2X2 + cixi^ 

£^8 xf + x| + C5X1X| + C4X1X| + C3X| + C2X| + C1X1X2 — S2X2 — SlXl 

(sx^ + C5X2 + C4X2 + C1X2 , 5x2 + 3C5X1X2 + 2C4X1X2 + 3c3x| + 2c2X2 + cixi^ 

Table 2 The left column indicates the A,D,E type of the Coxeter-Dynkin diagram or generic caustic. The right column 
lists the associated universal local forms of the smooth (n — l)-parameter family of general functions Fc,s, along with their 
(n — 3)-parameter family of induced general maps fc between planes — see the two-component expressions (... , ...). The 
given classification is due to Arnold 1973 Credits: Table from 



10 Universal Local Magnification Relations: Generic Caustics Beyond Codimension 3 

10.1 The Infinite Family of A, D, E Caustics 

In Arnold's classification of stable simple Lagrangian map-germs of n-dimensional Lagrangian subman- 
ifolds by their generating family Fc,s [6] (also, see Arnold, Gusein-Zade, and Varchenko [3 p. 330-31]), 
he found a deep connection between his classification and the Coxeter-Dynkin diagrams of the simple 
Lie algebras of types An (n > 2), D„ {n > 4), E^, E^, Eg. This classification is shown in Table [2] [3] and 
is known as the A, D, E classification of caustic singularities. 

The generic caustic singularities up to codimension 5 are given as follows in the Arnold A, D, E 
notation, where the numbers nD in parentheses indicate the codimension: 

1. (ID) A2 is a fold. 

2. (2D) A3 is a cusp. 

3. (3D) A^ is a swallowtail, an elliptic umbilic, and a hyperbolic umbilic. 

4. (4D) Ar, is a butterfly, D5 a parabolic umbilic. 

5. (5D) Aq is a wigwam, a 2nd elliptic umbilic, D'^ a 2nd hyperbolic umbilic, Eq symbolic umbilic. 

Remark: Up to codimension 5, all Lagrangian maps can be approximated by stable Lagrangian map- 
germs [6,. However, for codimension 6 or higher, this is no longer the case; unstable mappings form an 
open dense set [6ji8j. 

10.2 Universal Magnification Relations for the Family of A, D, E Caustics 

In 2009, Aazami and AOP [3] proved a univeral local magnification relation theorem for generic general 
mappings between planes exhibiting any caustic singularity appearing in Arnold's A, Z), E family: 

Theorem 16 [3] For any of the generic smooth (n — 1) -parameter family of general functions i^c,s {or 
induced general mappings fc) in the A, D, E classification, and for any non- caustic point s in the indicated 
region, the following results hold for the magnification 9Jti = 5iri(xi ; s) ; 

1. An {n > 2) obeys the magnification relation in the n-image region: X)r=i ~ ^' 

2. Dn (n > 4) obeys the magnification relation in the n-image region: Yl^=i ~ 

3. Eq obeys the magnification relation in the six-image region: ~ 0' 

4- E-j obeys the magnification relation in the seven-image region: X)I=i ~ 0, 
5. Eg obeys the magnification relation in the eight-image region: ~ 0- 
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Remark: Theorem does not follow directly from the Euler- Jacobi formula, the multi-dimensional 
residue integral approach [20], or the Lefschetz fixed point theory method [81], because some of the 
singularities have fixed points at infinity. 



10.3 On the Proof of Theorem[T6| 

The proof of Theorem [T6| given in [3] employed the Euler trace formula. This formula was shown by 
Aazami and AOP in 2009 to be a corollary of a more general result they established about polynomials: 

Theorem 17 [2] Let ip{x) = anX^^ + • • • + aix + oq G 'C[x\ be any polynomial with distinct roots Xi, and 
let h{x) ^ R he any rational function, where R is the suhring of rational functions that are defined at the 
roots of ip{x). Let 

/i*(a;) — Cn-ix"^^ + • • • + cix + co 
be the unique polynomial representative of the coset h{x) G R/{ip{x)) and let 

r{x) = bn-ix"^'^ H ^bix + bo 

be the unique polynomial representative of the coset ip'{x)h(x) € R/{ip{x)). Then the coefficients of r{x) 
are given in terms of the coefficients of h^{x) and (p{x) by the following recursive relation: 

bn-i = Cn-lbn-i,n-l + ' ' ' + Ci6„_i_i + Co6„_i,o i = 1, . . . , 71 , 

with 

' bn-ifl ^ (n- (i - 1)) a„_(i_i) , i = 1, . . . , n , 

bn-i,k = bn-i,k-l + ^'n-(4+l),fc-l : i — I, . . . ,n , k — 1, . . . ,n — 1 , 

fin 

where b-i^k-i = 0. 



Corollary 3 (Euler Trace Formula) Assume the hypotheses and notation of Theorem [7^ For any 
rational function h{x) G R, the following holds: 

" b 
V/i(2:i) = 
^ an 

See [10] for a residue calculus approach to the Euler trace formula. 
One can now show that the total signed magnification satisfies: 

y^, = (16) 

^ an 

For all of the caustic singularities appearing in the infinite family of A„ (n > 2), D„ {n > 4), EgjE-j, E% 
singularities, the coefficient was shown to be zero, and Theorem 1161 was thereby proved. We will 
illustrate the method of proof here in the case of the hyperbolic umbilic. See [2l[3] for a detailed treatment. 
The induced map fc corresponding to the hyperbolic umbilic is given by: 

^c{xi,X2) = (-32:^ - CX2 , -Sxl - CXi). 

Let s = (si, S2) be a target point lying in the four-image region. The four lensed images of s are obtained 
by solving for the equation 

{-Sxj - CX2 , -3x2 - cxi) = (si, 52)- (17) 

To use the Euler trace formula in the form (|16p . we begin by eliminating X2 to obtain a polynomial in 
the variable xi: 

/^{xi) = —Ssf — c?S2 — c^xi — 18six^ — 21x\. 



22 



The magnification of a lensed image of s under fc is 9}t(xi, X2) = l/(— + 36xia:2)- To convert tliis into 
a rational function in the single variable xi, we substitute for X2 via (|17p to obtain: 

— — 3dsiXi — 108xJ 

A direct calculation now yields: 

(p'{xi)m{xi) = c. 



Thus the unique polynomial representative in the coset ^p'{xl)^JfJl(xl) is the polynomial r{xi) = c (in the 
notation of Theorem [T71 VJl{xi) = h{xi)). Since 6„_i = 63 = 0, (|16p tells us immediately that 

9Jli + 9JI2 + M3 + 3714 = 0. 
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